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Abstract — I report on work done with Tomislav Prokopec (Heidelberg, 
CERN). We computed the one loop self-energy of a massless fermion dur- 
ing inflation. When the fermion is free it experiences only a small amount 
of inflationary particle production owing to the conformal invariance of its 
classical action. However, when the fermion is Yukawa coupled to a mass- 
less, minimally coupled scalar, there is copious production of fermions. In 
a more complicated model this effect might generate baryon asymmetry 
during inflation. 

1 Introduction 

Tomislav Prokopec and I have recently computed the one loop self-energy, 
during de Sitter inflation, of a massless fermion which is Yukawa coupled to 
a massless, minimally coupled scalar. The diagrammatic representation of 
what we did is given in Fig. ^ This would be a trivial exercise in flat space, 
and a fairly uninteresting one. The inflationary background geometry makes 
the calculation both nontrivial and physically interesting. What we found is 
that inflation engenders copious production of scalars which then decay into 
fermion-anti-fermion pairs. The first fact was known but the second is novel 
and potentially quite significant. In a more complicated theory it may provide 
a mechanism for generating a baryon asymmetry during inflation. 
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Figure 1: The one loop fermion self- energy. Fermion lines have an arrow, 
scalar lines do not. The final diagram gives the contribution of field strength 
renormalizat ion . 
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Most of the theoretical technology we used was developed in a previous 
computation, with Ola Tornkvist, of the one loop vacuum polariztion induced 
by massless scalar QED during inflation P3I2]- However, the result is very 
different. Whereas the photon develops a mass, which highly suppresses the 
production of photons during inflation, we find copious production of Yukawa 
coupled fermions. I will therefore devote most of this talk to explaining the 
physics in very simple terms. In Section 2 I give a brief review of inflation. The 
largest part of the talk is Section 3. In it I explain why the direct production 
of massless, minimally coupled scalars occurs during inflation, and why there 
is no significant direct production of fermions. Section 4 presents the result, 
without giving the derivation, and shows how it implies fermion production. 

2 Inflation 

On the largest scales the universe is amazingly homogeneous and isotropic 
0. It also seems to be devoid of spatial curvature |4j. The spacetime geometry 
consistent with these three features is characterized by the following simple 
invariant element, 

cis^ = —dt'^ + a?{t)dx ■ dx . (1) 

The coordinate t represents physical time, the same as it does in flat space. 
However, the physical distance between x and y is not given by their Euclidean 
norm, \\x — y\\, but rather by a(t)\\x — y\\. Because it converts coordinate 
distance into physical distance a{t) is known as the scale factor. 

Although the scale factor is not directly measurable, three simple observ- 
able quantities can be constructed from it, 

c^n ^ ^ d , , ad _ H 

a{tj a ii^ 

The redshift z gives the proportional increase in the wavelength of light emitted 
at time t and received at the current time, to. Redshift is often used to measure 
cosmological time, even for epochs from which we detect no radiation. The 
Hubble parameter H{t) gives the rate at which the universe is expanding. It's 
current value is, Hq = (71 ~ 2.3 x IQ-^^ Hz @]. The deceleration 
parameter q(t) is less well measured. Observations of Type la supernova are 
consistent with a current value of go — ~-6 [SJ. 

Inflation is defined as a phase of accelerated expansion, that is, q(t) < 
with H{t) > 0. From the current values of the cosmological parameters 
one can see that the universe seems to be in such a phase now. However, I 
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wish to discuss primordial inflation, which is conjectured to have occurred at 
something hke 10~^^ seconds after the beginning of the universe with a Hubble 
parameter 55 orders of magnitude larger than it is today. There are many 
reasons for believing that the very early universe underwent such a phase jH]. I 
will confine myself to reviewing how inflation resolves the smoothness problem. 
This can be summed up in the question, why does the large scale universe 
possess such a simple geometry (jT))? 

To understand the problem we need to compare the distance light can travel 
from the beginning of the universe to the time of some observable event, with 
the distance it can travel from then to the present. From the invariant element 
(jl)) we see that the radial position of a light ray obeys, dr = ±dt/a{t). The 
minus sign gives the past light- cone of the point = (to, 0), whereas the plus 
sign gives the future light- cone of a point = {ti, 0) at the beginning of the 
universe, 

J, _ dt _ fio,s dt 



'tabs a(t) Ju a{t) 

We can observe thermal radiation from the time of decoupling {zdec — 1089) 
whose temperature is isotropic to one part in 10^. Unless the universe sim- 
ply began this way — which seems unlikely — equilibrium must have been 
estabhshed by causal processes. In other words, we must have -Rfuture > -Rpast- 
Suppose that, during the period ti < t < t2, the deceleration parameter 
is constant q{t) = qi. In that case we can obtain explicit expressions for the 
Hubble parameter and the scale factor in terms of their values at t = ti, 

= TXTTT^VTTT a(t) = ai[l + . (4) 

i + [i-\-qi)ni(t-ti) L J 

These expressions permit us to evaluate the fundamental integral involved in 
the past and future light-cones ©, 

Jti a[t) QiHiqi'- J h qi 1 02/32 aiHi) 

Although go is negative, this is a recent event (2; ~ 1) which followed a 
long period of nearly perfect matter domination with g = +|. Much before 
the time of matter-radiation equality {zeq — 3200) the universe was almost 
perfectly radiation-dominated, which corresponds to g = +1. To simplify the 
computation we will ignore the the recent phase of acceleration and also the 
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transition periods, 

{Vl + 2; y Z < Zeq 

The cosmic microwave radiation was emitted within about a hundred redshifts 
of Zdec < Zeg, SO the past hght-cone is, 

_2fl 112 

The future hght-cone derives from both epochs and it depends shghtly upon 
the beginning redshift, z^eg, 

R^^^^^^ ^ J_ i ^ _ 1 1 + _±_ _ V^T+^ I (8) 

One maximizes -Rfuture by taking Zheg — 00, but it isn't enough. Under 
the assumption of g = +1 before z^g we are forced to conclude that the 
2- dimensional surface wc can see from the time of decoupling consists of 
(-Rpast /-Rfuture)^ — 2200 regious which cannot have exchanged even a photon 
since the beginning of time! So how did they reach equilibrium? 

This embarrassment resulted from the fact that the upper limit of integra- 
tion dominates i?future for positive deceleration. Inflation solves the problem 
by positing a very early epoch of negative deceleration. This makes the lower 
limit of -Rfuture dominate, so the future light-cone can be made as large as 
necessary by increasing Zbeg 



3 Inflationary Particle Production 
3.1 Virtual particles in flat space 

To understand why certain kinds of virtual particles can be ripped out 
of the vacuum during inflation it is instructive to review what constrains the 
lifetimes of virtual particles in flat space. A particle with wave vector k and 
mass m has energy, 

E(k) = ^Jm^ + \\k\\^ . (9) 

According to quantum field theory, virtual particles are continually emerg- 
ing from the vacuum with all different momenta. This process can conserve 
momentum if a virtual pair emerges with opposite 3-momenta, but it cannot 
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conserve energy. Before the pair came into existence the energy was zero, af- 
terwards it was 2E{k). The energy-time uncertainty principle asserts that we 
cannot detect this violation provided the pair annihilates in less than a time 
At given by the inequality, 

t+At ^ I 

dt'2E{k)<l =^ At<—-^ (10) 

2E[k) 

Of course the formalism of quantum field theory automatically incorporates 
these effects, without the need to invoke additional principles. However, the 
energy-time uncertainty principle allows one to understand many things. For 
example, virtual electron-positron pairs polarize the vacuum the most because 
they are the lightest charged particles and hence have the longest time to align 
themselves with applied fields. Similarly, it is long wave length virtual particles 
that survive the longest, which is one way of understanding long range forces. 

3.2 Virtual particles in cosmology 

In a homogeneous and isotropic geometry particles are still labeled 
by their constant, co moving wave vectors. However, because k involves an 
inverse wave length it must be divided by the scale factor in computing physical 
quantities such as the energy. 



E{t,k) = ^m^ +\\k\\ya^{t) . (11) 

Let us now reconsider the emergence of a virtual pair of such particles with 
wave vectors ±fc. If they emerge at time t, the energy-time uncertainty prin- 
ciple implies we will detect no violation of energy conservation provided they 
annihilate within a time At such that, 

t+At ^ 

dt'2E{t',k) <1 . (12) 



What can increase At? Obviously anything that reduces E{t', k). At fixed 
wave vector k and time t' this is accomplished by taking the mass to zero. Zero 
mass simplifies the integrand in (jl2j) to 2||/c||/a(t). Up to a constant this is the 
same expression as in the light-cones Q! We have already seen that negative 
deceleration (and hence inflation) causes these integrals to be dominated by 
their lower limits. Hence they cannot become arbitrary large as At goes to 
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infinity. The most negative deceleration consistent with stabihty is = — 1, 
for which the uncertainty bound gives, 



i-t+At ^ 9\\k\\ r n 

m = and a(t) = aic""^' =^ / dt' 2E{t, k) = 1 - e"^''^* . (13) 

Jt Hia[t) L J 

We therefore conclude that any virtual particle pair which emerges with < 
Hia{t) can survive for ever. 

3.3 Conformal invariance is the kiss of death 

One might imagine that the big obstacle to inflationary particle production 
is nonzero mass, but this is not so. At the enormous energy scales envisaged 
for primordial inflation all the known particles are effectively massless. The 
reason most of them still do not experience inflationary particle production is 
that they possess a symmetry suppressing the rate at which virtual particles 
emerge from the vacuum. 

The killer symmetry is known as conformal invariance. It rescales the fields 
which carry scalars (0), spin ^ fermions {ip), vectors (A^^) and gravitons {g^u) 
by powers of an arbitrary function of space and time fl{x), 

(j)(x) n-\x)(l){x) , A^{x) A^{x) , (14) 
ip{x) ^r^x)i!{x) , g^,uix) n\x)g^y{x) . (15) 

A simple, conformally invariant theory is electromagnetism, 

Cem = --^F^pFp^g^Pg^''^ , (16) 

where F^^ = d^A^ — d^A^. Since the vector potential is unaffected by a confor- 
mal transformation, the Lagrangian's invariance follows from the combination 
of inverse metrics and the square root of the determinant, 

g-Pg^^^ - n-'g'-P ■ Q-'g^^ ■ = g'^'g^'^V^g . (17) 

The connection between conformal invariance and cosmology derives from 
the existence of a coordinate system in which the general homogeneous and 
isotropic metric ((T)) becomes just g^u = o^Vt^y^ 



dt = adrj 



ds^ = a^i^—drj^ + dx ■ dxj = a^rj^j^^dx^dx^ . (18) 
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If we take a conformally invariant Lagrangian and express it in terms of the 
conformally rescaled fields with fl = 1/a, it is just as if the theory was in 
flat space! Some examples are the conformally coupled scalar, massless Dirac 
fermions, and electromagnetism, 

>Cccs = -ld^cl>d,cl>g^''V^ - ^cP'RV^ = -^d,{acl>)d,{acl>)n^'' , (19) 
/:Dirac = ?^eV0'9^-^^McdJ^')^y^ = ia^)Ytd^{ah) , (20) 



Physics does not depend upon local, invertible field redefinitions, so a con- 
formally invariant theory cannot be different, in conformal coordinates, than 
it is in flat space. This applies to all physical quantities, including the rate 
at which virtual particles emerge from the vacuum, per unit conformal time. 
Converting back to physical time gives the stated suppression, 

dn dri dn 1 / \ 

— = —— = - nat space rate . ( 22 ) 

dt dtdr] ^ ^ ^ 

Therefore any conformally invariant, massless virtual particles which happen 
to emerge with < Hia{t) will become real, but not many will emerge. 

3.4 Massless minimally coupled scalars 

Most familiar particles become (classically) conformally invariant when 
their masses are taken to zero. The exceptions are gravitons and massless, 
minimally coupled scalars. The Lagrangian density for the later is, 

>Cmmcs = ~d^<l>d,<f>g''''y^ = ~a''d^<f>d,<f>r'' . (23) 

Integrating over space to obtain the Lagrangian, then using Parseval's theorem 
to convert to Fourier space, gives a form we can recognize, 

d^xCuMCS^^J ^{a'^(t)|0(t,fc)|'-a(t)||^f |0(t,fc)P}. (24) 

— * 

Each wave vector k represents an independent harmonic oscillator with a time 
dependent mass, m{t) ~ a^{t), and frequency, u!{t) — \\k\\/a{ty. 
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For the case of de Sitter inflation {a{t) = aiC^'^) the mode functions take 
a simple form, 



Hi ik \ JK. 



(j){t, k) = u{t, k)a{k) + k)a)[—k) where u{t, k) = — ^===(l Jj)'' 

\ ^ (25) 

To understand the meaning of the operators a{k) and a"'" (A;), note that the 
minimum energy in wave vector k at time t is uj{t). However, the state with 
this energy does not evolve onto itself. Indeed, this theory has no stationary 
states! A reasonable "vacuum" is the state that was minimum energy in the 
distant past. This is known as Bunch- Davies vacuum and it is defined by 
a{k)\n) = 0. 

For Bunch-Davies vacuum the 0-point energy in wave vector k is, 

Eo{t,k) = la\t)\u{t,k)\^ + la{t)\m^^^^ (26) 
/ / za[i) 4||k|| 

The first term on the far right represents the irreducible, minimum energy. The 
second term gives the extra energy due to inflationary particle production. 
Since the energy of a particle of wave number k is ||/c||/a(t) we can easily 
compute the number of particles. 

As one might expect from the preceding discussion, this is much less than one 
at very early times, and it becomes order one when the wave number just 
begins to satisfy the uncertainty bound 

3.5 Scalar decay to fermions 

Although we have seen that inflation produces enormous numbers of mass- 
less, minimally coupled scalars, the conformal invariance of the Dirac La- 
grangian ^Id^ implies that there can be no comparable, direct production of 
fermions. However, it is still possible to make lots of fermions during infla- 
tion by allowing the scalars to decay into them. This can be accomplished by 
making scalars and fermions interact through a Yukawa coupling, 

^Yukawa = -^^MMCS + -^^Dirac — F^^/'t/^a/ — (28) 

= -la^d^(Pd,(Pr]^''+{Ji^)-f^'id^{ah)-Ta''fPij. (29) 
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Figure 2: Scalar decay into a fermion- anti- fermion pair. The initial state 
scalar has wave vector k. The final state fermions have wave vectors pi = p 
and p2 = k — p. 

The decay process is depicted in Fig. |21 Just as with infiationary particle 
production, 3-momentum is conserved but energy is generally not, 



k = p+ {k — p) 



AE{t) = \\p\\ + \\k-p\\ - \\k\ 



a{t) ■ 



(30) 



In fiat space the decay would be only be allowed for the infinitesimal phase 
space in which three momenta are co- linear. However, just as for particle 
production, the infiationary expansion weakens the constraint imposed by the 
energy-time uncertainty principle so as to permit the decay over a large volume 
of phase space. 



t+At 



dt'AE{t') < 1 =^ 
This is the physics behind our result. 



\p\\ + \\k-p\\ - \\k\\ <Ha{t) . (31) 



4 The Calculation 



I now specialize to the de Sitter scale factor, a{t) = e = —l/Hrj. It is 
useful to express the propagators in terms of the following conformal coordinate 
interval. 



Ax^fx; x') 



S'l|2 



(Ir^-VI-^e)^ (32) 

Owing to the conformal invariance of the free Lagrangian (j^Uj) . the fermion 
propagator is a trivial rescaling of the fiat space result. 



iSij{x; x') = {aa') ^^^' iQ^ 



1 



Atc^ Ax'^ix; x' 



(33) 
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where a = a{t) and a' = a{t') are the scale factors evaluated at the two points. 
The free scalar is not conformally invariant so its propagator has an additional 
term, 



iA{x; x') 



1 



477^ 1 Ax fx; x' 



2 



(34) 



The (pipiipj vertex and the fermion field strength renormalization are, 

— iTa'^6ij and i5Z2{aa')^'^ijid^5'^{x — x') . (35) 

Computing one loop diagrams in position space is easy. One simply mul- 
tiplies the various vertices and propagators. There are no integrations. The 
two diagrams of Fig. [T] give. 



— I 



x; X 



= {^—iV 5ikjiSM{x] x){^—iTa"^5(,j^iA{x] x') 

+t6Z2{aa')^^j id^5\x - x') . (36) 

Well, it isn't quite that easy! The actual calculation was done using the 
Schwinger-Keldysh formalism, which involves summing diagrams with slight 
variations in the ie prescription of expression \JjSi- 

One must also account for the ultraviolet divergences which manifest as 
products that are not integrable functions of x'^ and x''^. One deals with these 
by partially integrating to reduce negative powers of Ax^(x;x') until the re- 
sult is integrable. This procedure can be implemented so as to segregate the 
divergence to a delta function that can then be absorbed by 6Z2. Doing this 
rigorously requires the use of a regulator. We employed dimensional regular- 
ization, which makes slight changes in the fermion propagator ()33p and the 
vertices and major changes in the scalar propagator (jH^ Qj. 

Most of the formalism has been discussed in print ^2^ so I will confine myself 
to quoting the fully renormalized result. 



2%3 



2%2 



^ 7,^- idf.d'^ {e{Arj)e{A7]-Ax) (in [/i^ ( Ar/^- Ax^) 



ln(aa')(aa') 27^^*. id^6'^{x — x) 



d^d^ [e{Ar])e{AT]-Ax) \n[H\A7]^-Ax^)] } . (37) 



The first line is the conformally rescaled flat space result. Note that it involves 
a renormalization scale fi. The second line is the well-known contribution from 
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the conformal anomaly. The intrinsic de Sitter resuh is on the third hne. It is 
distinguished from the other two terms by its extra factor of aa'. 
The self-energy gives quantum corrections to the Dirac equation, 



a2ilji{x)) + I d^x' iEj {x\x')'4)j{x') = . 



(3^ 



Note that the flat space part of the one loop self-energy ()37|) has the same 
number of scale factors as the tree order differential operator. Perturbation 
theory only makes sense if the coupling constant is small, so we may assume 

^ 1. This means that the flat space part of the one loop self-energy can 
be ignored. So too can the conformal anomaly, which is only enhanced by a 
factor of ln(a). However, we cannot necessarily ignore the de Sitter correction. 
Although it is also down by F^, it is enhanced by a scale factor, which becomes 
enormously large during inflation. It therefore makes sense to keep only the 
de Sitter contribution of (jHTjl . 

Because the background is spatially translation invariant it also makes sense 
to look for plane wave solutions. 



ipi{ri,x) 



ik-x 



(39) 



Because the self-energy conserves helicity, we may as well also specialize to 
2-component helicity eigenstates, 

X^{V) = (^^1^] ) With k ■ aXL,R = ±kXL,R . (40) 
For the left-handed spinor the result of performing the spatial integrations is. 



= {zdo±k)xL{v) 

J^dv'a'xL{v')e^"'^'' 



21n(2iJAr/) + l+ / dr 

Jo T 



.(41) 



To obtain the result for a right-handed spinor with ± helicity, simply change 
the signs of ±fc, ^ikAt] and ±ir. 

Recovering the full time evolution of this sort of nonlocal mode solution re- 
quires numerical integration. However, it is straightforward to get the asymp- 
totic behavior for late times. In this regime the physical 3-momentum has 
redshifted to essentially zero, so we can simplify ()41|) by setting k = 0. This 
obliterates the distinctions between left and right handedness, and between 
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different helicities. The equation can be further simphfied by transforming 
from conformal time rj to co-moving time t, 



dtX + 



dt'x{ln 



. 



(42) 



Now factor the leading exponential out of the logarithm, 



In 



2e 



-Ht' 



2e 



-Ht 



-Ht' + ln(2) + In 



-HAt 



(43) 



It is an excellent approximation to retain only the term —Ht'. Acting another 
derivative results in a local equation, 

whose approximate solution can be obtained by the WKB method, 

^TH 



X-^iHt)-^exp —(Ht) 
L on 



(45) 
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